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Quark Masses on the Lattice: Light and Heavy 

V. Lubicz a 

a Dipartimento di Fisica, Universita di Roma Tre and INFN, Sezione di Roma III, 
Via della Vasca Navale 84, 1-00146 Rome, Italy 

I review the current status of lattice calculations of light and heavy quark masses. Significant progresses, in 
these studies, have been allowed by the introduction of improved actions and non-perturbative renormalization 
techniques. Current determinations of light quark masses are accurate at the level of 20%, where the main source 
of uncertainty is represented by the quenching error. The determination of the bottom quark mass is accurate 
at the impressive level of 2%. As final averages of lattice results, I quote m u[ j(2 GeV) = (4.5 ± 1.0) MeV, 
m s (2 GeV) = (110 ± 25) MeV and m b (m b ) = (4.26 ± 0.09) GeV. 
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1. INTRODUCTION 

Calculations of quark masses are one of the 
most intensive subject of investigation for lat- 
tice QCD. An accurate determination of these 
parameters is in fact extremely important, for 
both phcnomcnological and theoretical applica- 
tions. The charm and bottom masses, for in- 
stance, enter through the heavy quark expansion, 
the theoretical expressions of several cross sec- 
tions and decay rates. From a theoretical point of 
view, an accurate determination of quark masses 
may give insights on the physics of flavour, re- 
vealing relations between masses and mixing an- 
gles, or specific textures in the quark mass matrix, 
which may originate from still uncovered flavour 
symmetries. 

The values of quark masses cannot be directly 
measured in the experiments because quarks are 
confined inside the hadrons. On the other hand, 
quark masses are fundamental (free) parameters 
of the theory and, as such, they cannot be com- 
puted on the basis of purely theoretical conside- 
rations. The values of quark masses can be only 
determined by comparing the theoretical evalua- 
tion of a given physical quantity, which depends 
on quark masses, with the corresponding experi- 
mental value. Typically the pion, kaon and <p me- 
son masses are used on the lattice to compute the 
values of the light quark masses, whereas the b- 
quark mass is determined by computing the mass 



of the B or the T mesons. Different choices are 
all equivalent in principle, and the differences in 
the results, obtained by using different hadron 
masses as input parameters, give an estimate of 
the systematic error. 

As all other parameters of the Standard Model 
lagrangian, quark masses can be defined as ef- 
fective couplings, which are both renormaliza- 
tion scheme and scale dependent. A scheme 
commonly adopted for quark masses is the MS 
scheme, with a renormalization scale chosen in 
the short-distance region in order to make this 
quantity accessible to perturbative calculations. 
It is a common practice to quote the values of the 
light quark masses at the renormalization scale 
fi = 2 GeV, whereas the heavy quark masses are 
usually given at the scale of the quark mass itself, 
e.g. rrib{rnb)- This convention will be followed in 
the present review. 

At Lattice 99, no plenary talk was dedicated 
to review lattice calculations of quark masses. 
For this reason, I will mainly report on re- 
sults presented in the last two years. With re- 
spect to previous determinations, significant pro- 
gresses in controlling the systematic errors have 
been allowed by the extensive implementation 
of non-perturbative 0(a)-improvement and 
non-perturbative renormalization techniques j^, 
||. The effect of these tools, on lattice calcula- 
tions of quark masses, will be discussed in sec. B. 
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The recent determinations of light quark masses 
will be presented in sec. [|, and of the 6-quark 
mass in sec. |]. Some comments on the deter- 
mination of the charm quark mass will be done 
in sec. ||, where I will also summarize the lattice 
"world averages" of quark masses and present my 
conclusions. 

2. DEFINITIONS AND THEORETICAL 
PROGRESSES 

2.1. Quark Masses and 0(a)-improvement 

Non perturbative definitions of quark masses 
are provided by the chiral Ward identities (WI) of 
QCD H| . These definitions allow us to express the 
renormalized quark mass, in a given scheme and 
at a given renormalization scale, in terms of lat- 
tice renormalization constants and bare quanti- 
ties (action parameters, matrix elements or corre- 
lation functions). Three independent definitions 
have been proposed so far which are all equiva- 
lent in principle, in the sense that they lead to 
the same value of the renormalized mass. Let us 
discuss these definitions in some details. 

1) Vector Ward Identity definition: 

In terms of renormalized quantities, the chiral 
vector WI reads: 

VM%\P) = (mi(p) - ma(A«)) (a\S(nM W 

where \a) and \(3) are arbitrary external states. 
Vector current conservation, within the lattice 
regularization, implies a lattice version of the vec- 
tor WI, from which the following definition of the 
renormalized quark mass, in terms of the (bare) 
Wilson parameter, can be derived: 

TO (M) = Z s 1 (a[i)m = Zs 1 ^)-^ (j£ - jf-) (2) 

An 0(a) -improved definition is obtained with the 
simple replacement M: 



m — > m (1 + ab m m) 



(b m = b s /2) 



(3) 



2) Axial Ward Identity definition: 

In the Wilson formulation of lattice QCD, axial 
symmetry is explicitly broken. The axial WI can 
be then imposed only for renormalized quantities, 



and, in the case of degenerate masses, takes the 
form: 



V ll (a\A ll \0) = 2m(n) (a|P( M )|/?) 



(4) 



From this identity, an alternative definition of the 
quark mass is derived: 



m(n) = 



Za to awi = — Yli^\At 
Zp(a/i) Zp{ajjL) 2 (a\P\i 



(5) 



In this case, C(a)-improvement is achieved with 
the replacement 



rn 



AWI 



[1 + a (b A - bp) m] ■ 
V^<q|(A^ + oc^V M P)[ 
2{a\P{n)\(3) 



(6) 



In the absence of explicit chiral symmetry brea- 
king induced by the Wilson term in the action, 
the vector and axial vector WI definitions are not 
independent anymore. This is the case, for in- 
stance, of staggered fermions and of the proposed 
realizations of Ginsparg- Wilson fermions (domain 
wall and overlap fermions). 

3) Definition from the Quark Propagator: 

This relatively new method has been proposed 
by APE [[| . It is based on the study of the quark 
propagator, computed in a fixed gauge, at large 
momentum. In this limit, an OPE can be derived 
for the trace of the quark propagator which, up 
to known logarithmic corrections, reads Q : 



Tr S{p) 



3p 4 



(W) + 0(1//) 



(7) 



The leading term of this expansion is proportional 
to the renormalized quark mass. The Wilson co- 
efficient of this term is fixed by an axial WI which, 
in RI-MOM scheme, implies B: 



£- x (p,p) 



(8) 



Note that the relevant renormalization constant, 
in this case, is the quark field renormalization 
constant Z q . Cancellation of leading 0(a) ef- 
fects is obtained by improving the quark field. 
Since this field is a non gauge-invariant operator, 
and the quark propagator is an off-shell correla- 
tion function, improvement involves mixing with 
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both non gauge-invariant operators and opera- 
tors vanishing by the equation of motion. The 
O(a)-improvement of the quark field is obtained 
through the replacement j7|: 

q — » (1 + ab q m) [1 + ac! q {Jj) + m ) + ac NG i^% (9) 



one finds, in the large momentum limit and up to 
small logarithmic corrections || : 
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Figure 1 . Large momentum behaviour of the lat- 
tice quark propagator. The upper figure shows 
the large O(a) term which has to be subtracted 
before extracting the subleading contribution pro- 
portional to the quark mass. The behaviour of 
the inverse quark propagator, after the subtrac- 
tion, is shown below. The results have been ob- 
tained by APE j^j using the non-perturbatively 
O(a)-improved Wilson action at (3 = 6.2. 

In the determination of quark masses with the 
propagator method, the improvement procedure 
is a crucial requirement. The reason is that the 
OPE of the quark propagator on the lattice is 
dominated by a large lattice artifact. From eq. (||) 



12^ S{P) 



-2a c' + 2ac NG iZ q 



a c a 



(10) 



The lattice quark propagator does not vanish at 
large p 2 , as dictated in the continuum limit by 
eq. (0), but it is dominated by terms of 0{a). 
This is illustrated in fig. |l|. The numerical study 
of the propagator at large p 2 , in the chiral limit, 
provides the value of the coefficient c q , which can 
thus be subtracted before computing the sublead- 
ing contribution to the OPE proportional to the 
quark mass. 

The improved calculations of quark masses, 
with the vector and axial-vector WI methods, re- 
quire the knowledge of several coefficients, ca, 
bA, bs and bp, besides the coefficient csw of 
the clover term in the action. Non-perturbative 
determinations of these coefficients have been 
performed in refs. P]|p§- One finds that the 
effect of improvement, in the determination of 
quark masses, is quite significant, even for light 
quarks. In the case of the strange quark mass, 
for instance, an estimate of the 0(a) effects, 
with unimproved Wilson fermions, has been per- 
formed by Bhattacharya and Gupta jl0| at Lat- 
tice 97, based on the analysis of the continuum 
extrapolation. They find that discretization er- 
rors, at two typical values of the lattice cou- 
pling, (3 = 6.0 and [3 = 6.2, are of the order 
of 25% and 18% respectively. These estimates 
can be compared with those obtained with non- 
perturbatively C(a)-improved Wilson fermions, 
from a recent study by ALPHA-UKQCD gl|. 
In this case, residual discretization effects are re- 
duced at the level of 7% and 3% respectively. Im- 
plementation of C(a)-improvement is therefore a 
crucial ingredient to significantly increase the ac- 
curacy of lattice calculations of quark masses. 

2.2. Non-perturbative Renormalization 

Another potential source of systematic errors, 
in lattice calculations of quark masses, is intro- 
duced by the truncation of the perturbative ex- 
pansion in the evaluation of the quark mass renor- 
malization constant. Depending on the defini- 
tion considered for the bare mass, this constant 
is given by Zg 1 , with the vector WI method, 
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Table 1 

Values of the strange quark mass, rn s (2 GeV), in MeV, as obtained by computing the quark mass 
renormalization constant either in one loop (boosted) perturbation theory (PT) or using the FJ-MOM 
non-perturbative method (NP). The relative difference between the two determinations is also shown. 



Action 





m s [PT] 


m s [NP] 


A 




NP-CLOVER 


6.2 


95(10)t 
81(9)* 


109(11) 
111(13) 


15% 
37% 


APE 98 fn§ 


KOGUT-SUSSKIND 


6.0 
6.2 
6.4 


80(4) 
81(6) 
83(4) 


114(5) 
110(8) 
108(5) 


43% 
36% 
30% 


JLQCD 99 @ 


DOMAIN WALL 


6.0 


106(8)(14)* 


130(11)(18) 


23% 


RBC 99 pi 



t From vector WI; * from axial WI. * The error is my estimate from the RBC data. 



Z^/Zp, with the axial WI method, or the quark 
field renormalization constant Z q with the prop- 
agator method. In perturbation theory, all these 
quantities are known only at one loop. In the 
last two years, however, in most of lattice quark 
masses calculations, the systematic error associ- 
ated with the evaluation of the renormalization 
constant has been reduced to a negligible amount 
by the use of non-perturbative renormalization 
techniques 

Non-perturbative renormalization is reviewed 
by Sint at this conference Jl^] . For this reason, I 
will only discuss here those particular aspects of 
this topic which are more closely related to the 
determination of quark masses. 

To show the impact of non-perturbative renor- 
malization in lattice calculations of quark masses, 
I present in table |l| the values of the strange 
quark mass, as obtained by computing the quark 
mass renormalization constant either in one-loop 
(boosted) perturbation theory or with the RI- 
MOM non-perturbative method of ref. Q. For 
illustrative purposes, the results of three calcula- 
tions are considered, which employed three diffe- 
rent fermionic actions: non-perturbative Clover, 
Kogut-Susskind and domain wall fermions. As 
can be seen from table |l|, the difference between 
the results obtained with perturbative and non- 
perturbative renormalization is rather large. De- 
pending on the action, the value of the coupling 
constant and the specific renormalization con- 
stant (Zs or Zp for Wilson fermions), the re- 
lative error varies in the range between ~ 15% 



and ~ 40%. Moreover, at least in the case of 
staggered fermions, this error decreases as one 
get closer to the continuum limit, but at a very 
slow rate. Therefore, non-perturbative renorma- 
lization is required to improve the accuracy of lat- 
tice calculations of quark masses at a significant 
level (below 15%). 

I believe that a more detailed discussion de- 
serves the non-perturbative calculation of the 
pseudoscalar renormalization constant, Zp, in the 
RI-MOM scheme. Although rather technical, this 
discussion is of relevance for the determination of 
quark masses with the axial WI method. Within 
the non-perturbative approach of ref. 0], Zp is 
evaluated from the amputated correlation func- 
tion of the pseudoscalar density, computed be- 
tween external off-shell quark states of momen- 
tum p. This function, Tp(p), once conveniently 
projected onto the matrix 75, is related to the 
trace of the inverse quark propagator by the ax- 
ial WI: mYp{p) = TrS , - 1 (p). The OPE of the 
quark propagator, given in eq. (0), thus implies: 



Tp(p)^C 1 (p 2 ) + C 2 (p 2 ) 



mp 2 



+ 0(l/ P 4 ) (11) 



The renormalization constant Zp is determined 
from the coefficient function C\. Eq. ( jll"| ) shows 
that the first power correction, which vanishes at 
large momentum, is however divergent in the chi- 
ral limit. Physically, this divergence is introduced 
by the coupling of the pseudoscalar density to the 
pion field, which is the Goldstone boson of the 
spontaneously broken chiral symmetry. It has 
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been emphasized, particularly in ref. |lq| , that, 
because of the singular behaviour of this term 
in the chiral limit, its contribution to the cor- 
relation function Tp(p) may not be suppressed 
enough at large p 2 , at least for those values of p 2 
in actual calculations at which discretization ef- 
fects are still under control. This may affect the 
determination of Zp and, according to ref. Jl6| , 
this systematic error may in fact be the reason of 
the large discrepancy between the perturbative 
and the non-perturbative determination of Zp. 
Since Z m ~ Zp 1 , the observation of ref. Q is of 
relevance for lattice calculations of quark masses. 

I believe that the warning of ref. |ll| deserves 
some consideration and a careful investigation. 
I would like to mention, therefore, three differ- 
ent methods which can be used to verify whether 
the contribution of the Goldstone pole, in the ac- 
tual calculation of Zp, is suppressed at the re- 
quired level. To study this problem, I will use 
the same APE data which have been used for the 
non-perturbative calculation of the renormaliza- 
tion constant of bilinear quark operators (Zy, Za, 
Zs, Zp and Zt) with the non-perturbative Clover 
action, at f3 = 6.2 |ft7| . The first method is based 
on the study of the scale dependence of the ratio 
Zp/Zs, computed non-perturbatively in the RI- 
MOM scheme. In the large momentum region, 
where power corrections are expected to be negli- 
gible, this ratio should exhibit a plateau, since the 
two renormalization constants have equal anoma- 
lous dimensions. This ratio is shown in fig. |[ 
and in the largc-p 2 region a plateau is clearly vi- 
sible. This suggests that all power corrections, 
including the contribution of the Goldstone pole, 
are suppressed enough. The second method is 
based on the study of the renormalization group 
invariant combination Zp(a/j,)/C(afj,), where the 
running factor C(a/j), computed in perturbation 
theory at the NLO, expresses the predicted lo- 
garithmic scale dependence of the renormaliza- 
tion constant. Again, in the region where non- 
perturbative power corrections are suppressed 
enough, the ratio Zp(a/j,)/C(afi) should behave 
as a constant. Indeed, this constant behaviour is 
observed in fig. ^. Finally, the most efficient way 
to get rid of the Goldstone boson contribution in 
the calculation of Zp (but the method can be ea- 
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Figure 2. The ratios Zp/Zg, Zp(aix)jC{a[i) and 
Zp(a[j,)/C(afj,) (see text for details) as a function 
of the scale (a/i) 2 . 

sily adapted to the calculation of other renormal- 
ization constants) has been recently suggested in 
ref. [Q. The idea is to consider a new definition 
of Zp, which I will denote with Zp, involving two 
flavours of non-degenerate quarks: 



Mm) 



m\ — 7712 



(12) 



Z q {n) miT P (fi, mi) - m 2 Tp(fj,, m 2 ) 

On the r.h.s. of eq. (|T^), the limit 777 1 ,m 2 — ► 
is understood. For 777 2 = (and at large values 
of the scale), the above definition reduces to the 
standard PJ-MOM definition of Zp. However, 
one immediately realizes that, at the leading or- 
der in the quark mass, the Goldstone boson con- 
tribution of eq. ([ll]) cancels in the definition of 
Zp, which is thus wellj^ehaved in the chiral limit. 
To the calculation of Zp, the warning of ref. JTg[ 
does not apply anymore. The ratio Zp(fj,)/C(afi) 
is plotted in fig. as a function of the scale. One 
can see, in this case, that the desired constant be- 
haviour is reached much faster than with the ratio 
Zp(fi) /C(afi). Moreover, the two ratios converge 
to the same value at large scales, within less than 
5%, an uncertainty which is of the same order of 
the other systematic errors involved in the cal- 
culation. This analysis supports the assumption 
of an adequate suppression of the Goldstone pole 
contribution in previous calculations of Zp. 

A significant accuracy in the non-perturbative 
determination of the quark mass renormalization 
constant has been achieved by ALPHA 0] , within 
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the Schrodinger functional approach. Besides de- 
termining the value of mass renormalization con- 
stant at a low hadronic scale, ALPHA has also 
computed, non-perturbatively, the running factor 
C(/x) over a large region extending from the low 
hadronic scale up to a very high scale. This cal- 
culation, which is based on a recursive finite-size 
scaling technique, allows a fully non-perturbative 
determination of the renormalization group in- 
variant quark mass. The non-perturbative calcu- 
lation of C{fj) turns out to be in very good agree- 
ment with the predictions of NLO perturbation 
theory. The advantage of this approach is that 
the matching factor, relating the renormalization 
group invariant mass to the MS mass, is known 
at N 4 LO in perturbation theory Jl|,|§. The re- 
maining perturbative uncertainty in this determi- 
nation is thus reduced at a completely negligible 
level. The ALPHA calculation of the mass renor- 
malization constant has been discussed by Suit 
at this conference, and I refer to his review for 
further details 

3. LIGHT QUARK MASSES 

3.1. Quenched Calculations 

The determination of light quark masses is, at 
present, one of the main fields of activity of lat- 
tice QCD simulations. This is illustrated in fig. ^, 
which presents a compilation of all lattice results 
for the strange quark mass, in the quenched ap- 
proximation, starting from the first lattice predic- 
tion for m s obtained with a NLO accuracy pi) . 
These results have been produced by using dif- 
ferent actions, different orders of improvement, 
different renormalization techniques and different 
choices of the input parameters, used to fix both 
the scale and the strange quark mass itself (ei- 
ther mx or m^). In spite of these differences, we 
find good consistency among the several determi- 
nations. In addition, all the results, with the pos- 
sible exception of the rather large value obtained 
by SES AM |||] , suggest an interval of allowed va- 
lues for the strange quark mass which is smaller 
than the one quoted by the particle data group 
(PDG) Q, indicated by the band in the figure. 
One can also observe that the lattice determina- 
tions of the strange quark mass cluster close to 
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Figure 3. Results of all lattice calculations of the 
strange quark mass obtained, in the quenched ap- 
proximation, since 1994- The band shows the es- 
timate of m s quoted by the particle data group. 
The dashed line indicates the unitarity lower 
bound on rn s . 

the unitarity lower bound, fn s ^ 100 MeV p3], 
obtained by using unitarity and the positivity of 
the pseudoscalar spectral function.^] Finally, the 
results which are plotted on the same vertical line 
in the figure have been obtained by using either 
the kaon or the <f> meson mass as input parameter. 
It has been shown that a systematic difference of 
approximately 20% is found between the two de- 
terminations within the quenched approximation. 

In order to derive an average value for the 
strange quark mass, I will concentrate in the fol- 
lowing on the more recent and accurate results, 
which are given in table || and shown in fig.^j with 
full circles. In this case, only the results obtained 
by using rriK as experimental input have been 
considered. 

The comparison among the results presented 
in table || becomes more significant if one tries 
to correct them by taking into account for the 
different systematic errors involved in the calcu- 

1 The N 2 LO bound of ref. slightly decreases to m s <; 
95 MeV, when the N 3 LO perturbative corrections to the 
pseudoscalar spectral function are taken into account |G5| . 
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Table 2 

Quenched lattice results for the strange quark mass obtained by using 
mass. Details concerning the action, the value of the lattice spacing, the 
the perturbative accuracy are also given. 



Action 



APE 98 |13| NP-Clover 
JLQCD 99 fll KS 
CP-PACS 99 |6| Wilson 
ALPHA-UKQCD 99 NP-Clover 

QCDSF 99 @ NP-Clover 

APE 99 (| NP-Clover 

CP-PACS 00 H) MF-Clover 

RBC 00 M DWF 



in input the value of the kaon 
rcnormalization procedure and 

PT m s {2 GeV)/ McV 



N 2 LO 
N 2 LO 
NLO 
N 4 LO 
N 4 LO 
N 3 LO 
NLO 



7 GeV 
-> 
-> 
-> 
-> 
7 GeV 
-> 



K) 

(m p ) 
(m p ) 
{ro,fic) 
(ro) 

(»»p) 



NP-RI 
NP-RI 

PT 
NP-SF 
NP-SF 
NP-RI 

PT 



1.9 GeV (m p ) NP-RI N 3 LO 



111 ± 12 
106 ±7 
116 ±3 
97 ±4 
105 ±4 
111 ±9 

noil 

110 ±2 ± 22 t 
105 ± 6 ± 21 * 



From vector WI; * from axial WI. 



QUENCHED 



PDG 

APE 98 
JLQCD 99 

CP-PACS 99 

ALPHA + 
UKQCD 99 

QCDSF 99 
APE 99 
CP-PACS C 
RBC 00 
AVERAGE 
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rio 130 

m B (2 GeV) 



150 
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Figure 4. Values of the strange quark mass 
obtained from recent lattice calculations in the 
quenched approximation (full circles). Empty cir- 
cles denote the values and errors obtained by the 
author attempting to correct for the different sys- 
tematics. 

lations. In this way, I obtain the estimates shown 
in fig.|] as empty circles. 

All results in table || have been obtained adopt- 
ing a non-perturbative renormalization techni- 
que, with the only exceptions of the two CP- 
PACS determinations |2^,^| , in which the quark 
mass renormalization constant has been evalu- 



ated using one-loop perturbation theory. An ad- 
ditional uncertainty, due to the use of perturba- 
tion theory, should then be added to these results, 
to account for the corresponding systematic er- 
ror. This uncertainty, however, is difficult to be 
quantified, because two different renormalization 
constants [Zg 1 and Za/Zp) have been used in 
the calculation and, moreover, non-perturbative 
determinations of renormalization constants with 
the Iwasaki action used in ref. |p8| have not been 
performed yet. In fig. |], a systematic error of 10% 
has been added to the CP-PACS results. This er- 
ror, however, may be underestimated, according 
to the analysis of the previous section (see ta- 
ble @). 

The values of the strange quark mass by AL- 
PHA-UKQCD and QCDSF @ have been 
obtained by fixing the lattice spacing from the 
value of the phenomenological parameter ro |$2| , 
rather than from m p as in all other determina- 
tions. QCDSF finds that the choice of ro is 
roughly equivalent to using m p , while ALPHA- 
UKQCD finds that this is equivalent to setting 
the scale from f K . ALPHA-UKQCD also esti- 
mates that a 10% higher value of rn s would have 
been obtained by fixing the scale from the value 
of the nucleon mass. Since the scale determined 
from vrip turns out to be typically in between the 
values obtained from fx and rnjv, in order to 
compare with the other determinations I have in- 
creased by 5% in fig.[| the values of m s obtained 
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by ALPHA-UKQCD and QCDSF. 

The other corrections to the results presented 
in table p| ar e smaller. The APE determinations 
of refs. ||5[|lS|1. obtained with the non-perturbative 
Clover action at (3 = 6.2, have not been extrapo- 
lated to the continuum limit. In this case, the 
more extensive analysis by ALPHA-UKQCD Jll[] 
shows that the value of the strange quark mass is 
underestimated by approximately 3%. Finally, in 
the calculations by APE Jlf] and JLQCD the 
conversion from the non-perturbative RI-MOM 
renormalization scheme to the MS scheme has 
been done at the N 2 LO in perturbation the- 
ory |30[ , since at the time when these studies 
have been performed the N 3 LO perturbative cal- 
culation of |nj was not available yet. In fig-H in 
order to account for the difference between N LO 
and N 3 LO, the results of refs. Q and Q have 
been decreased by 3%. 

I would like to emphasize that the previous dis- 
cussion, besides exploiting the different sources 
and sizes of systematic errors, also shows that a 
rather good level of statistical and systematic ac- 
curacy has been achieved in the determinations of 
light quark masses within the quenched approxi- 
mation. Indeed, the values of m s shown in fig. || 
are in remarkable agreement within each other, 
and even more when the differences in the sys- 
tematics have been taken into account. From the 
spread of the results, and accounting for the typi- 
cal uncertainty on the scale in the quenched case 
(~ 10%), I obtain as a final average of the lattice 
results, within the quenched approximation, 

to, (2 GeV) QUEN = (110 ±15) MeV (13) 

This value is shown in fig|4|, together with the se- 
veral lattice results and the average value of to s 
quoted by the PDG |^3| . Note that the uncer- 
tainty in eq. (|l^) is approximately three times 
smaller than the one quoted by the PDG, al- 
though, in the former, the effect of the quenching 
approximation has not been taken into account 
yet. 

In the light quark sector, the other quantity, 
besides to s , which is accessible to lattice calcu- 
lations is the average value of the up and down 
quark masses, m u d = (m u + rrid) /2. Rather then 
the value of to U£ j itself, however, I find more con- 



Table 3 

Quenched lattice results for the ratio of the 
strange to the average up-down quark masses. 
2m s /(m u + m d ) 

APE 98 |J 24.6 ±2.2 

JLQCD 99 pi 25.1 ± 1.7 

CP-PACS 9 9M 25.3 ± 1.0 

QCDSF 99 jfl 23.9 ±1.1 

APE 99 U 23.4 ±2.4 

CP-PACS 00 H 25.2 ±1.0 
The errors are my estimates based on the original 
data. 

venient to consider the ratio of the strange to 
the average up-down quark masses, because in 
this ratio many of the systematic uncertainties, 
coming for instance from the renormalization con- 
stants or the choice of the scale, are expected to 
partially cancel. The values of m s /m u d from the 
most recent lattice calculations are collected in 
table and shown in fig. H. The quoted errors 
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• — 


— 1 


CP-PACS 99 


1 — • 


H 


QCDSF 99 


— • — 1 




APE 99 | 






CP-PACS 00 


• 1 




ChPT 











15 20 25 30 

2m s /(m u +m d ) 



Figure 5. Ratio of the strange to the average up- 
down quark masses as obtained from recent lattice 
calculations in the quenched approximation (cir- 
cles). The errors are my estimates based on the 
original data. The lowest point (diamond) repre- 
sents the prediction of chiral perturbation theory. 

are my estimates based on the original data, be- 
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cause the value of this ratio has not been quoted 
directly in the original papers. 

The relevant result, shown in fig. [|, is that the 
lattice predictions for m s /m u d are in very good 
agreement with the value 24.4 ± 1.5 predicted by 
chiral perturbation theory [f33| . By using this in- 
formation and eq. (|l^) , I obtain the quenched ave- 
rage: 

m ud {2 GeV) QUEN = (4.5 ± 0.6) MeV . (14) 

Finally, I would like to mention that the first 
(quenched) lattice calculation of m u d with over- 
lap fermions has been presented at this confer- 
ence [Q. The preliminary result is m u d{2 GeV) 
=5.5(5) MeV, in good agreement with the avera- 
ge ©. 

3.2. Unquenched Results 

Several unquenched calculations of light quark 
masses, with two flavours of dynamical quarks, 
have been presented at this conference. A compi- 
lation of the results, together with details of the 
simulations, is given in table || With the excep- 
tion of the SESAM f| and CP-PACS @ deter- 
minations, all other results presented in the table 
are still preliminary. 

A striking result shown in table ^ is the value 
of the ratio of the strange to the average up- 
down quark masses obtained by SESAM | p5[ : 
m s /m u d ~ 55. This value is larger, by roughly a 
factor 2, than all other quenched and unquenched 
lattice determinations (including the quenched 
result by SESAM [p2[), and in disagreement with 
the prediction of chiral perturbation theory. The 
reason for such a large value of m s /m u d obtained 
by SESAM is the way in which the unquenched 
analysis has been performed. Indeed, as discussed 
below, by following a different procedure they ob- 
tain m s /m u d — 30 [^2|, which is much closer to 
the expected value. I believe this point is worth 
to be discussed in details, because it is also of 
relevance for other unquenched studies. 

In numerical simulations with dynamical fer- 
mions all physical quantities, as well as the lattice 
spacing and the critical value of the hopping pa- 
rameter, are functions of the two parameters en- 
tering in the action, namely the coupling constant 
and the sea quark mass. Theoretical considera- 



tions suggest that the dependence of the lattice 
spacing on the value of the sea quark mass should 
be exponential ]Icj| . Indeed, the lattice scale de- 
pends exponentially on the effective coupling con- 
stant which, in turn, is affected by the values of 
quark masses because quarks run in the loops. A 
large dependence of the lattice spacing on the sea 
quark mass is in fact observed in actual calcula- 
tions. For instance, JLQCD (3^] finds that, in 
the range of K sea used in the simulation (corre- 
sponding to mps/mv — 0.6 — 0.8), the inverse 
lattice spacing varies approximately by 30%, be- 
tween ~ 1.5 and ~ 2.0 GeV. The important point 
is that such a dependence, if not properly taken 
into account, may lead to uncontrolled systematic 
effects in the various extrapolations, in both the 
sea and valence quark masses. 

In order to control these effects, an appropriate 
procedure for the unquenched analysis has been 
advocated [^,^l). It consists in performing all 
the analysis, including the calibration of the lat- 
tice spacing and the calculation of the relevant 
physical quantities, at each fixed value of K sea . 
The simulations performed at fixed K sea can be 
thought of as pseudo- quenched simulations, which 
come closer to the description of the real world as 
the sea quark mass approaches its physical value. 
Only when all the quantities, computed at fixed 
K sea , are expressed in physical units, the results 
can be extrapolated in the sea quark mass. 

This procedure has been followed by APE |37j 
in the calculation of the strange quark mass. 
At two fixed masses of the sea quark, both the 
strange and the average up-down quark masses 
have been computed. The results are shown 
in fig. ^, in terms of the ratios m u d/mK* and 
rUs/niK*- Presenting these ratios is equivalent 
to show the quark masses in physical units. The 
values of m u d are then extrapolated in the sea 
quark mass, also computed in physical units, up 
to the point corresponding to m sea = m u d- Once 
the physical value of m u d has been determined in 
this way, a similar extrapolation is performed for 
the strange mass. From fig. |^ one can see that 
the final extrapolations in the sea quark mass are 
indeed smooth.n 



2 Since, however, only two values of K aea have been simu- 
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Table 4 

Simulation details and physical results of unquenched lattice calculations of light quark masses. 





Action 


a" 1 






m s (2 GeV) 


m s /m ud 


r~\ T TT71 "NT i T TM/^» 

m QUEN^UNQ 


SESAM 98 H 


Wilson 


2.3 GeV 


4 


PT 


151(30) 


{ m K,<f>) 


55(12) 


1.10(24) 


MILC 99 || 


Fatlink 


1.9 GeV 


1 


PT 


113(11) 
125(9) 


{m K ) 
(m ) 


22(4) 


1.08(13) 


APE 00 Hi 


Wilson 


2.6 GeV 


2 


NP-PJ 


112(15) 
108(26) 


{m K ) 
( m <p) 


26(2) 


1.09(20) 


CP-PACS 00 |f 


| MF-Clover 


a 


12 


PT 




(m K ) 
i m 4>) 


26(2) 


1.25(7) 


JLQCD 00 |§ 


NP-Clover 


2.0 GeV 


5 


PT 


94(2)t 
88(3)* 
109(4)* 
102(6)* 


{m K ) 
( m <p) 






QCDSF + 
UKQCD 00 (H 


NP-Clover 


2.0 GeV 


6 


PT 


90(5) 


{m K ) 


26(2) 





* From vector WI; * from axial WI. The errors on the ratios m s /m u d and m® 13 ™ /m™Q are my estimates 
based on the original data. 



In the analysis by SESAM and CP-PACS the 
pseudoscalar and vector meson masses have been 
fitted simultaneously in both the sea and the va- 
lence quark masses, with all quantities expressed 
in lattice units. The large dependence of the lat- 
tice spacing on the sea quark mass may not be 
taken properly into account in this way. Indeed, 
SESAM finds that the value of m s /m u d changes 
by a factor 2 when the analysis, instead, is per- 
formed at fixed values of K sea , as advocated be- 
fore. Moreover, the result obtained with the anal- 
ysis performed at fixed K sea is much closer to the 
expected value of this ratio. The systematic ef- 
fects, introduced by the neglecting of the lattice 
spacing dependence on the sea quark mass, are 
difficult to evaluate a priori, and depend on the 
specific values of the parameters used in the cal- 
culation. JLQCD, for instance, attempted the 
chiral extrapolations by using either dimension- 
less or dimensionful quantities |38| (although con- 
sidering different orders in the expansions), and 



latcd by APE, the results should be taken as preliminary. 



obtained for m s compatible results. 

Since unquenched studies of light quark masses 
have not yet reached the same degree of accu- 
racy achieved in quenched calculations, in order 
to obtain an estimate of the quenching effect it is 
convenient to consider directly the ratio between 
quenched and unquenched results. Some of the 
systematic errors are expected to cancel in this 
ratio, when both determinations are performed 
by using the same action, the same renormaliza- 
tion procedure and the same choice of physical 
inputs. The results for the ratio m% VEN /M™ Q 
are presented in table ^. In the very exten- 
sive calculation performed by CP-PACS the value 
TO QUEN / / m UNQ = 1.25(8) is obtained, suggesting 
a sizable decrease of the quark mass in the un- 
quenched case. This result, however, is not con- 
firmed by the other (less extensive) studies by 
SESAM, APE and MILC, which find a decrease, 
in the unquenched case, rather of the order of 
10%. Given these results, and bearing in mind the 
uncertainties in the analysis methods, I believe 
that, in order to quote a final estimate of light 
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0.05 7 0.70 
m sea (3 GeV)/m K , 



Figure 6. Extrapolation of the average up- down 
(a) and the strange (b) quark masses, computed at 
fixed K sea , to the physical value of the sea quark 
mass, m sea = m u d- 



quark masses, it is more appropriate to include 
the quenching error as an additional systematic 
uncertainty in the quenched averages (eqs. (13) 
and ([u])), rather than varying the central values. 
Assuming this error to be of the order of 20 MeV 
in the case of the strange mass, I then obtain: 



m s (2 GeV) = (110 ± 15 ± 20) MeV (15) 
m ud (2 GeV) = (4.5 ± 0.6 ± 0.8) MeV (16) 

Note also that in the unquenched case (see ta- 
ble U) the differences between the determinations 
of the strange quark mass from ttik and m^, if 
any, are rather small. 

4. THE b-QUARK MASS 

The o-quark cannot be directly simulated on 
the lattice since its mass is larger than typical va- 
lues of the ultraviolet cutoff in present lattice cal- 
culations (a -1 ~ 3 GeV). The 6 mass, however, is 
also larger than the typical energy scale of strong 
interactions. Thus, the heavy degrees of freedom 
of the 6-quark can be integrated out. Beauty 
hadrons may be therefore simulated on the lattice 
in the framework of a low-energy effective theory. 
Indeed, in the past years, many lattice simula- 



tions of Heavy Quark Effective Theory (HQET) 
and Non-Relativistic QCD (NRQCD) have been 
performed. 

Within the effective theory, the mass of a B- 
meson, Mb, is related to the pole mass of the b- 
quark. Up to higher-order 1/mb corrections, this 
relation has the form: 



pole 



5m 



(17) 



where e is the so-called binding energy, which 
can be computed non-perturbatively by a nume- 
rical simulation, and dm is the residual mass, ge- 
nerated by radiative corrections, which can be 
evaluated in perturbation theory. The perturba- 
tive calculation of the residual mass and the non- 
perturbative (lattice) calculation of the binding 
energy thus allow a determination of the b pole 
mass. The result can be then translated into the 
MS mass, nib, by using perturbation theory. 

An important observation, concerning this pro- 
cedure, is that the binding energy e is not a phy- 
sical quantity, and it is affected by power diver- 
gences proportional to the inverse lattice spacing, 
1/a. These divergences are canceled by similar 
singularities in the residual mass 8m. Moreover, 
the pole mass m^ ole , and consequently 8m, are 
also affected by rcnormalon singularities, which 
introduce in their definitions an uncertainty of 
the order of i\.qcd ^^f^f • These singularities are 
then canceled by the perturbative series relating 
the pole mass and the MS mass [^4|. As a result, 
the MS mass, nib, is a finite, well defined, short- 
distance quantity. In the actual calculation, how- 
ever, since the residual mass is computed up to 
a finite order in perturbation theory, only a par- 
tial cancellation of both power divergences and 
the renormalon singularities occurs. For this rea- 
son, it is crucial to compute the residual mass 8m 
up to the highest possible order in perturbation 
theory. Moreover, because of the incomplete can- 
cellation of power divergences, the lattice spacing 
in such calculations cannot be taken too small. 

At present, the most accurate determination of 
8m has been obtained in the framework of HQET. 
In terms of the lattice bare coupling constant, the 
perturbative expansion can be written as: 



a 8m = C\ a + C 2 + C 3 al + 



(18) 
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The two- loop coefficient, C2 has been derived 
analytically, by Martinelli and Sachrajda 45 , 



from the study of the asymptotic behaviour of 
Wilson loops. In the quenched case (Nj — 0), 
also the coefficient C3 has been evaluated, by 
Di Renzo and Scorzato [Q, using the method 
of numerical stochastic perturbation theory. An 
important check of this approach is obtained 
from the comparison of their numerical results 
for the first two coefficients, C\ = 2.09(4) and 
C2 = 10.7(7), with the analytically known values: 
C x = 2.1173 and C 2 = 11.152 JIB). The value ob- 
tained for the third coefficient, C3 = 86.2(5), is 
also consistent with the determination by Lepage 
et al. [J47j|, C3 = 81(2), obtained, within a com- 
pletely different approach, by fitting the results 
of small coupling Monte Carlo calculations. 

These combined theoretical efforts allow a de- 
termination of the 6-quark mass which is accu- 
rate, in the quenched approximation, up to the 
N 3 LO. Two independent results have been ob- 
tained so far: 

_ QUEN = ( 4 30 ± Q Q5 ± aQ5 j QeV j 48 j £ 19 j 

m b QUEN = (4.34 ± 0.03 ± 0.06) GeV @ (20) 

nicely in agreement within each other. Remar- 
kably, the two determinations have been derived 
using different approaches. The APE result ( |l9| ) 
has been obtained within HQET, while the deter- 
mination j20|), by Collins et al, is obtained from 
the NRQCD study of the B-meson spectrum in 
the static limit. The last error in eqs. ( |l9| ) and 
( p(i| ) represents the residual uncertainty due to the 
neglecting of higher orders in the perturbative ex- 
pansion of 5m. For comparison, this uncertainty 
was estimated to be of the order of 200 MeV at 
NLO §§ and 100 MeV at N 2 LO |§ respectively. 

Results compatible with those in eqs. ( |l9| ) and 
( |20| ) have been obtained from the study of the B 
and the T spectrum with NRQCD [|l],||. In this 
case, however, only the coefficient C± is known, 
so that the achieved accuracy is at NLO. These 
results are thus affected by a larger perturbative 
uncertainty, which can be estimated, on the basis 
of a renormalon analysis [Q, to be of the order 
of 100 MeV. 

The first unquenched calculation of the 6-quark 
mass has been performed by APE this year |40|. 



The perturbative accuracy is at N 2 LO, since the 
three-loop coefficient, in the expansion of dm, is 
yet unknown in the unquenched case. Their re- 
sult, 



m b = (4.26 ± 0.06 ± 0.07) GeV . 



(21) 



represents to date the most accurate determina- 
tion of the &-quark mass from lattice QCD calcu- 
lations. Remarkably, the relative uncertainty is 
reduced at the level of 2% (90 MeV). 

A preliminary analysis ]49| ] performed in the 
static limit with NRQCD indicates a decrease of 
r?Th, in the unquenched case, by approximately 70 
MeV. Combined with the quenched determina- 
tion (|2(]) , this result leads to an estimate of the 
6-quark mass in good agreement with (^l|). 

5. CONCLUSIONS 

Significant progresses in lattice calculations of 
light and heavy quark masses have been allowed 
by the introduction of improved actions and non- 
perturbative renormalization techniques. 

Current determinations of light quark masses 
have reached a good level of statistical and sys- 
tematic accuracy, which is of the order of 10% 
within the quenched approximation. The main 
source of uncertainty, in this case, is due to the 
quenched approximation. Several unquenched re- 
sults, with two flavours of dynamical quarks, have 
been presented at this conference. The effect of 
quenching is found to be in the range of 10-20%, 
but further investigations are required to obtain 
a firmer estimate. 

I have not reported, in this review, lattice cal- 
culations of the charm quark mass, since recent 
results for this quantity are missing. The charm 
mass is of great phenomenological interest since 
it enters, through the heavy quark expansion, 
the theoretical predictions of several cross sec- 
tions and decay rates. For this reasons, new lat- 
tice calculations with the non-perturbatively im- 
proved action and non-perturbative renormaliza- 
tion techniques is highly recommended. 

Enormous progresses have been achieved in lat- 
tice calculations of the bottom quark mass. The 
present accuracy is at N 3 LO in the quenched ap- 
proximation, and at N 2 LO in the unquenched 
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case. The relative uncertainty is reduced at the 
impressive level of 2%. 

Final averages of lattice calculations for quark 
masses have been given in eqs. (|ll|), (|l6|) and (pl|). 
By combining statistical and systematic errors, 
these correspond to: 

m ud {2 GeV) = (4.5 ± 1.0) MeV 

m s (2 GeV) = (110 ± 25) MeV (22) 

m b (m b ) = (4.26 ± 0.09) GeV . 
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